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ABSTRACT Density functional theory (DFT) and coupled cluster theory (CCSD(T)) calculations were 
employed to investigate the geometric and electronic structures of a range of dinuclear molybdenum sulfide 
clusters, Mo2S,” and MooS, (n = 4~~8). The results showed that the sulfur atoms tended to occupy the terminal 
sites of the clusters continuously in the process of sequential sulfidation. After the oxidation state of Mo atoms 
reached the maximum of +6, diverse disulfur ligands emerged in the sulfur-rich Mo2S,,’° (n = 7, 8) clusters. 
The driving forces of removing a sulfur atom from different S ligands in Mo2S,- (n = 4~8) clusters, 
especially from those disulfur units, were evaluated. The corresponding order may provide insight into the 
pretreatment of fresh MoS catalysts. Vertical detachment energies (VDEs) were predicted according to the 
Generalized Koopmans’ theorem, and then the photoelectron spectra (PES) were simulated. Molecular orbital 
and spin density values were analyzed to elucidate the chemical bonding and the evolutionary behavior in the 
dinuclear molybdenum sulfide clusters. 
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1 INTRODUCTION 


Molybdenum sulfide has intensely caught people's eyes due to their wide applications in the fields of 
mechanical treatments, catalysts, energy materials and so forth!'*!. In the case of catalysts, it is known that 
molybdenum sulfide can serve as a promising low-cost alternative to platinum or other noble metals catalyst in 
hydrogen evolution reaction (HER)?""4!, Additionally, MoS>-based catalysts have been widely used in the 
process of petroleum refining!!, including the hydrodesulfurization (HDS) and hydrodenitrogenation (HDN) 
reactions. Meanwhile, considerable attention has been paid to study the active sites in these catalysts!!°*5!. It is 
normally accepted that the sulfur vacancies (or coordination unsaturated sites (CUS)) of MoS2 are the key 
active centers!”° 771. However, Besenbacher et al.'*7! also put forward that the specific brim sites without sulfur 
vacancies are also responsible for the HDS activity. Furthermore, Vrubel et al.'''! showed that the reduced 
molybdenum sulfides containing the disulfide ligands appeared to be catalytically active. Several molybdenum 
complexes which contain various S2 ligands have been synthesized recently, which were considered to be able 
to mimic the MoS> edge sites for the catalytic hydrogen generation'**3"!_ Despite the progressive work, more 
detailed studies on the active sites of molybdenum sulfide catalysts are still required. 

Previously, many studies have been focused on the synthesis and characterization of a variety of 
molybdenum polysulfido complexes and their precursors in the condensed phase, and various disulfide species 
were obtained in these complexes'?**° 3351. Among them, the [Mo2S7]* anion was found to possess a 
terminal S2% ligand!” The [Mo2Ss]*~ was considered to own two terminal S2% ligands!??!. The [Mo2S 12]*> was 
reported to contain four terminal S2% and two bridging S2% ligands°!. The [Mo3S13]* was predicted to have 
three bridging S2% and three terminal S2% ligands!”*!. These molybdenum complexes containing the persulfido 
(S2*-) ligands may lead to the corresponding complexes with supersulfido (S27) ligands via futher electron 
transferation, but less is known about their supersulfido complexes°~*"!. As is known to all, gas-phase clusters 
can serve as effective molecular models gaining microcosmic understanding of the sophisticated surface 
structures and the catalytic processes at the molecular level], Theoretical calculations have played an 
indispensable role in describing the accurate structure and properties of gas-phase clusters!'“*!, Over the past 


few years, considerable efforts have been devoted to studying the mono- and multi-nuclear molybdenum 


sulfide gas-phase clusters'7-!, Infrared spectra combined with DFT calculations of neutral mono-nuclear MS, 
M = Cr, Mo, W; n = 1~3) clusters have been studied by Andrews's group”. Joint experimental 
photoelectron spectra and theoretical investigations on a variety of mono- and multi-nuclear molybdenum 
sulfide clusters have also been reported by Gemming et al.?- 57 581, Jiao et al. presented the structure and 
reactivity of Mo3So cluster which can be taken as a model for amorphous molybdenum sulfide MoS3'*4!. An ab 
initio study on the structural stability of Mo-S clusters and the size specific stoichiometries of magic clusters 
was reported by Murugan et al.°!, wherein Mo2Ss consisting of two bridging S and three terminal S atoms was 
expected to be the magic cluster. In spite of these efforts, the systematical theoretical investigations on the 
gas-phase molybdenum sulfide clusters are still required. 

In our previous work'*!: ©, we have reported a theoretical study on the mono-nuclear molybdenum 
sulfide clusters, MoS, “° (n = 1~6). To mimic the geometric and electronic properties of molybdenum sulfide 
surfaces and defects, larger Mo,S, clusters in size may be interesting. In the present work, extensive density 
functional theory (DFT) and coupled cluster theory (CCSD(T)) calculations were performed to elucidate the 
structural and electronic properties of a range of di-nuclear molybdenum sulfide clusters, Mo2S,, and Mo2S, (n 
= 4-8). The current study represents our continuous research interest in various clusters aiming at providing 
well-defined molecular models for bulk surfaces and catalysts'®!°*!, According to our calculations, a behavior 
of structural evolution was found with the exception of neutral Mo2Ss. The neutral Mo2Sg can be viewed as 
replacing two of the terminal S atoms in Mo2S¢ by the same number of S2 units. Interestingly, diverse disulfur 
ligands, including the supersulfido (S27) ligands, emerged in sulfur-rich clusters Mo2S,,° (n = 7, 8). It was 
found that the disulfur species may have a key impact on the catalytic activity’ !. Our calculations showed 
that the reduced reactions for removing one sulfur atom from the sulfur-rich clusters were spontaneous, 
whereas the reactions were nonspontaneous if the sulfur atom was removed from the sulfur-deficient clusters. 
The results suggested that the S2 units may play an important role in removing the sulfur atoms from the edge 


sites of fresh MoS: catalysts. 


2 COMPUTATIONAL METHODS 


The calculation details for this study are similar to our earlier studies on mono-nuclear molybdenum 
sulfide clusters, MoS, ”° (n = 1~6)'°! 61, Density functional theory (DFT) calculations employing the B3LYP 
hybrid functional %7? were carried out using Gaussian 03 program!’?!. B3LYP was widely used in quantum 
chemistry!™4!, Furthermore, B3LYP had been applied in the other Mo-S systems, which showed good 
agreement with the experimental data>!:>4.°°!, Additionally, B3LYP also gave reasonably good results which 
were compared to the available experimental data in our previous work'*! ©, As discussed below, we used the 
results with B3LYP functional for further discussion. A host of initial structures considering different spin 
states and geometric symmetry were evaluated, and the search for the most stable structures was first 
performed using the triple-¢ valence plus polarization (def2-TZVP) basis set!’>77! and the corresponding 
Stuttgart effective core potential for Mo?! (denoted as L-BS hereafter). Then the selected low-lying isomers 
(AE < 0.50 eV) were further re-optimized at the B3LYP level with the larger basis sets, i.e., the Stuttgart 
relativistic small core basis set and efficient core potential!’® 7! augmented with two f-type and one g-type 
polarization functions (C(f) = 0.338, 1.223; &g) = 0.744) for molybdenum"®! and the aug-cc-pvtz basis set for 
sulfur and hydrogen'*!*?! (denoted as H-BS hereafter). Scalar relativistic effects were taken into account via 
the quasi-relativistic pseudo-potentials. Vibrational frequency was calculated at the same level of theory to 
confirm that the reported minima have no imaginary frequency. The relative stabilities of several energetically 
close-lying isomers (AE < 0.40 eV) were further distinguished with the help of higher-level CCSD(T)!****! 
single point calculations with the H-BS basis sets at the B3LYP optimized geometries. Vertical electron 
detachment energies (VDEs) were calculated on the basis of the generalized Koopmans’ theorem!®*! which had 
been described detailedly in our previous studies'®!®!, All DFT calculations were performed using the 
Gaussian 03 software package. The CCSD(T) calculations were performed with the MOLPRO 2010.1 


package!"!, The frontier molecular orbitals were visualized using the VMD software!!! 


3 THEORETICAL RESULTS 


The optimized geometries for the ground-state and selected energetically low-lying isomers of MoS, 


and MooS, (n = 4~8) at the B3LYP/H-BS level of theory are displayed in Figs. 1~5. Their relative energies 
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including those isomers within 0.40 eV at the B3LYP/H-BS level together with the results of single-point 
CCSD(T) calculations are collected in Table 1. Alternative optimized results at the B3LYP/L-BS level for 
MoS,” (n = 4~8) are given in the Supporting Information (Figs. S1~S5). 
3.1 Sulfur-deficient clusters: Mo2S, and MoS, (n = 4, 5) 

Previously, the tribridged structure with three bridging S atoms and the dibridged structure with two 


bridging S atoms and one terminal S atom were reported to be the possible ground states of the Mo2S3 
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cluste To search for the ground states of Mo2S4“, a host of initial structures were taken into 
consideration, including the above mentioned dibridged and tribridged strucutres. Based on the calculations, a 
nonplanar dibridged structure (C;, 7A") is shown to be the lowest-energy structure of Mo2S4 (Fig. 1a). Another 
triplet state 7A2) with higher symmetry C2, is only 0.04 eV (Fig. 1b) higher in energy. The corresponding 
quintet state (B1) with C2, symmetry is located to be 0.15 eV (Fig. 1c) higher in energy. The quintet state (B) 
with Cr, symmetry is 0.21 eV (Fig. 1d) higher in energy. The difference between these two quintet states is 
that two terminal S atoms in the former show a syn relationship, whereas the latter has a anti configuration. In 
addition, the triplet state (Cs, °A") with the anti configuration is 0.30 eV (Fig. le) higher in energy. It seems 
that the syn configuration is more stable than the anti configuration. The previous studies'4”: 48 56. 581 have also 
supported the dibridged structure which has the syn configuration to be the ground state of Mo2S4. For the 
anionic species, the dibridged structure (C2, B1) with syn configuration is found to be the lowest-energy 
structure of Mo2S4 (Fig. 1g). Another two dibridged isomers, C2; (4B) and C2, (7B1), are located to be 0.18 eV 
(Fig. 1h) and 0.34 eV (Fig. 1i) higher in energy, respectively. Other optimized isomers for both the neutral and 
the anion are much higher in energy and thus not listed in the current paper. More optimized structures at the 
B3LYP/L-BS level are available in the Supporting Information (Fig. S1). 

Starting from the ground states of Mo2S4~” clusters, extensive structural searches revealed the singlet (Cs, 
IA’) to be the ground state of Mo2Ss (Fig. 2a). It can be viewed as adding a terminal S atom to the dibridged 
Mo2S4. The previous studies by Murugan et al.5®% and Gemming et al.*! also supported our results. The 
corresponding triplet state (Cs, °A’) is 0.17 eV (Fig. 2b) higher in energy. For the anionic Mo2S57, the ground 
state (Fig. 2c) is found to be a doublet state (Cs, 7A'), whose geometry is similar to that of the neutral ground 


state. A quartet (Ca, *A2) is located to be 0.22 eV (Fig. 2d) above the ground state. 
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3.2 Stoichiometric clusters: Mo2S¢ and MoS 

The ground state of Mo2S6 is identified to have C2, ('A1) symmetry (Fig. 3a). It can be deemed as adding a 
terminal S atom to the Mo2Ss ground state. Each Mo atom in MooS¢ is tetra-coordinated with two terminal S 
atoms and two bridging S atoms. This result is in consistent with the previous study by Gemming et al.!°!. It 
should be mentioned that the isomer with a terminal S2 unit was expected to be the most stable structure of 
Moo2S¢ by Murugan et al.°!. But this isomer is less stable (> 0.50 eV) than the structure in Fig. 3a according to 
our calculation (Fig. S3). For the anionic Mo2S67, an open-shell (Az) structure with D2, symmetry is found to 


be the ground state (Fig. 3b). 
3.3 Sulfur-rich clusters: MoS, and MoS, (n = 7, 8) 

To our knowledge, no other computational studies have been reported for the gas-phase Mo2S7-” and 
Mo2S3~ clusters. On the basis of our calculations on Mo2S¢”, various initial configurations with different 
spin multiplicities were studied in search of the ground states of Mo2S7. As shown in Fig. 4a, the 
lowest-energy structure of Mo2S7 is predicated to be closed-shell (1A) with C2 symmetry wherein a bridging S2 
unit appears. It can be regarded as replacing one of the bridging S atoms in Mo2S¢ by a bridging S2 unit. The 
S-S bond length of the S2 unit in Mo2S7 (Fig. 4a) is 2.086 A.A triplet state (Cs, 3A") with a terminal S2 unit 
(Fig. 4b) is located 0.26 eV above the ground state. However, this structure becomes the ground state of the 
Mo2S7" anion (Fig. 4d), whereas the structure with a bridging S2 unit (Fig. 4e) is 0.39 eV higher in energy. The 
Mo2S7" ground state can be seen as replacing one of the terminal S atoms in Mo2S¢ by a terminal S2 unit. The 
terminal S2 ligand in Mo2S7 (Fig. 4d) is attached in a side-on fashion, and the corresponding S-S bond length 
is 2.097 A. 

In the case of Mo2Ss-”, we found several close-lying isomers of Mo2Ss~”” near the lowest-energy structure. 
As shown in Fig. 5a, the ground state of Mo2Ss is found to be *B2 state with C2, symmetry, which can be 
viewed as replacing two terminal S atoms in Mo2S¢ by two terminal S2 units in a syn configuration. The S—S 
bond lengths in the S2 units are calculated to be 1.999 A, and the Mo-S> bond lengths are 2.468 A. Another 
triplet state (C2, *B,) which also has two terminal S2 ligands but in a anti configuration is 0.20 eV (Fig. 5b) 
above the ground state. In addition, the isomer which can be seen as replacing the remaining bridging S in 
Mo2S7 by a bridging S2 unit is much higher in energy (> 0.5 eV; Fig. S5). For the anion, the most stable 


structure of Mo2Sg” (Fig. 5c) is a doublet state with C2, symmetry, in which two terminal S2 ligands are in a 


syn relationship as in the case of the neutral. The anti isomer (C2, 7Bz) is shown to be 0.06 eV (Fig. 5d) higher 
in energy. Additionally, the isomer which contains a terminal S; ligand is 0.11 eV (Fig. 5e) higher in energy. 
Similar with the results of Mo2S4%, the calculations showed that the syn configuration seems to be more stable 
than the anti configuration for Mo2Sg™”. 
3.4 CCSD(T) single-point calculations for low-lying structures 

The low-lying isomers of Mo2S, (n = 4~8) and their anion species (within 0.40 eV at the B3LYP/H-BS 
level) were further evaluated using higher level single point CCSD(T) calculations at the B3LYP geometries. 
The relative energies of single point CCSD(T) calculations are summarized in Table 1. As a whole, the results 
of single point CCSD(T) calculations are in good line with that of DFIT/B3LYP calculations, except for the 
neutral Mo2S4. According to the CCSD(T) calculations (Table 1), the second lowest-energy isomer (Czy, 7A2; 
Fig. 1b) of Mo2S4 seems more stable than the isomer (Cs, 7A") shown in Fig. 1a. Although these two isomers 


are close in erengy (0.07 eV) at the CCSD(T) level, they are likewise similar in their geometries and electronic 


structures. Herein, the Mo2S4 (Czy, 7A2) shown in Fig. 1b is tentatively considered to be the ground state. 


4 DISSCUSION 


4.1 Interpretation of the simulated spectra and molecular orbital analyses 

The experimental PES spectra which can be used as an electronic “fingerprint” of a given cluster could 
provide valuable electronic information. On the basis of generalized Koopmans’ theorem, the vertical 
detachment energies (VDEs) for the identified Mo2S,” (n = 4~8) anionic ground states and the selected 
low-lying isomers (AE < 0.40 eV) were calculated (Table 2). The introduction to the principles of PES spectra 
simulation may be found in the early references!*”!, This simulation method has been extensively used in a 
number of previous studies'*!: © 94961 and show good agreement with the experimental spectra. The simulated 
PES spectra are presented in Fig. 12. The first vertical detachment energy (VDEis) trend as a function of S 
content (n) in MoS,” (n = 4~8) is displayed in Fig. 13. 
4.1.1 Mo2S4 and Mo2Sq- 


The frontier orbitals of Mo2S4 and Mo2S¢ are illustrated in Fig. 6. As mentioned above, the Mo2S4 (Cav, 
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3A) as shown in Fig. 1b is considered to be the ground state using the CCSD(T) single-point calculations. Its 
valence electronic configuration is (18a1)?(10b1)?(6a2)'(19a1)'. Addition of an electron into the empty 
antibonding orbital 13b2 of the neutral (Fig. 6a) would lead to the anionic ground state Mo2S4- (Cz, +B1) as 
shown in Fig. 1g. The corresponding valence electronic configuration is (10b1)?(18a1)?(6a2)!(13b2)'(19a1)!. 
Accordingly, the Mo—Mo distances increase from 2.616 to 2.760 A, as shown in Fig. 1b and 1g. The VDEis 
upon photodetachment from the singly occupied bonding orbital (19aı)!, which is mainly characterized by a 
Mo 4d orbital, is predicted to be 3.39 eV. The fully occupied 10b; MO and below are mainly of S 3p character. 
Other calculated VDEs from 1352 and below are presented in Table 2. 
4.1.2 MoSs and MoS 

The ground state of neutral Mo2Ss (Cs, !A') is predicated to be closed-shell species (Fig. 2a). The valence 
electronic configuration of Mo2Ss is (36a')?(37a')*. Addition of an electron into the empty bonding orbital 38a’ 
of the neutral (Fig. 7a) would lead to the anionic ground state Mo2Ss (Cs, ?A'; Fig. 2c) with a valence 
electronic configuration of (36a')?(37a')*(38a')!. As shown in Fig. 7b, the singly occupied molecular orbital 
(SOMO) 38a' and the fully occupied orbital 36a' primarily correspond to the Mo 4d orbitals, and the orbital 
37a’ is characterized by S 3p feature. The first detachment channel for Mo2Ss" is derived from the removal of 
the SOMO 38a’, for which the calculated VDE is 3.76 eV. 
4.1.3 MoS. and MoS 

The neutral Mo2S¢ (Fig. 3a) is stoichiometric in which each Mo achieves its highest oxidation state Mo. 
In other words, all of six valence electrons of Mo [4d°5s'] are used to form bonds with the S atoms. The 
valence electronic configuration of the neutral Mo2S6 (Ca, 'A1) cluster is (24a1)?(12b1)°(19b2)*. The highest 
occupied molecular orbital (HOMO) 19b2 and lowest unoccupied molecular orbital (LUMO) 25a; are depicted 
in Fig. 8a. All MOs from HOMO and below are S 3p-based orbitals, and the LUMO is primarily featured by 
the Mo 4d orbital. When one extra electron is added to the neutral Mo2S6, the ground state of Mo2Se6 (Fig. 3b) 
would still maintain the skeleton of neutral (Fig. 3a) but the rhombus Mo2S2 unit in the neutral is distorted 
relative to the anion. The valence electronic configuration of Mo2S6~ (Dan, 7Ag) is (4b3g)*(10b2u)?(15a¢)!. In this 
Dən structure, the extra electron delocalized on the two Mo atoms, as is evidenced from its SOMO 15a; (Fig. 
8b). The photodetachment from SOMO 15a, yields the VDEis with a calculated value of 4.51 eV. The 
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remaining calculated VDEs from 102, and below are presented in Table 2. 
4.1.4 MoS; and MoS; 

For the sulfur-rich species Mo2S7, the structure with a S2 unit in the bridging fashion is predicated to be 
the neutral ground state (Fig. 4a), for which the valence electronic configuration is (35a)?(33b)?(36a)>(34b). 
All the MOs are featured by 3p orbital of S atom (Fig. 9a). The ground state of anion Mo2S7 (Ci, 7A) is 
predicted to be open-shell with a terminal S2 group (Fig. 4d). The valence electronic configuration of this anion 
is (65a)*(66a)*(67a)*(68a)'(69a)*(70a)*(71a)*. The calculated S-S bond length in this S2 group is 2.097 A, 
which is similar with the free S2% ('X,*) dianion (2.180 A calculated at the same level). The two fully occupied 
orbitals 71a and 65a correspond to the z* orbitals of S2 group (Fig. 9b). Hence, the Mo2S7 Cı ('A) can be 
considered as the addition of a S2% unit to the cationic Mo2Ss5*. Photodetachment from the fully occupied 
orbital 71a yields the first PES band with the calculated VDEs of 4.32 eV (a) and 4.26 eV (£). The detachment 
origining from the singly occupid orbital 68a; requires a higher energy (5.00 eV). 

4.1.5 MoS; and Mo2Ss 

The frontier MOs of Mo2Sg and MooSs" are illustrated in Fig. 10. The ground state of Mo2Ss is found to 
be *B> state with C2, symmetry (Fig. 5a). The valence electronic configuration for the neutral species is 
(16b1)?(26a1)*(17b1)!(13a2)'(18b1)°(27a1)*. The corresponding frontier MO pictures are shown in Fig. 10a, in 
which two singly occupied orbitals 13a2, 17b, and two doubly occupied orbitals 16; and 18; correspond to 
the z* orbitals of two bound S2 units in Mo2Ss (C2, 3B2). Furthermore, the S-S bond length of the S2 moiety 
(1.999 Ain Fig. 5a) is very close to that of free S2 (Mg) anion (2.029 A calculated at the same level). The spin 
density analyses further certify that two unpaired electrons are separately located on two S2 units (Fig. 11a). 
Therefore, the Mo2Sg (Cz, *B2) may be viewed as two So units adhered to the cationic Mo»S,4**. For the 
anionic species, the valence electronic configuration is (16b1)?(26a1)?(17b1)?(27a1)°(13a2)°(18b1)!. As shown 
in Fig. 10b, the singly occupied MO 18b; and three doubly occupied MOs 13a2, 17b; and 16b; correspond to 
the z* orbitals of the S-S moieties. The bond length of S2 unit (2.055 A) is between that of the bound S2 and 
S2% anions (1.999 A in Fig. 5a and 2.097 A in Fig. 4d, respectively). The spin density analyses show that an 
unpaired electron is equiprobably shared by two terminal S2 groups (Fig. 11b). Therefore, the Mo2Ss7 (Cav, 7B1; 
Fig. 5c) cluster can be described as a resonance hybrid of two equivalent C, structures with both S27 and S2% 
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units. As shown in Table 2 and Fig. 12e, photodetachment from the fully occupied 13a2 orbital (Fig. 10b) of 
MooSs° yields the first PES band with the calculated VDEis of 3.58 eV (). 
4.2 Structural evolution of Mo2S,~” (n = 4~~8) clusters 

As noted above, a series of thiomolybdate dianion Mo2S,?> (n = 6~9) have been prepared and 
characterized in condensed phase!**!, It was pointed out that six homologues [Mo2S,]*> (n = 6~12) dianions 
may be obtained from any other by either adding sulfur or removing sulfur by triphenylphosphine (Ph3P)"°?!. In 
our paper, similar structural evolutions are found for the gas-phase Mo2S,,° (n = 4~8) clusters. 

It is worth to mention that Mo2S6 showed the highest stability among a series of Mo2S, clusters by 
Gemming et al.!58!, They indicated that the structure of Mo2S« may be used as a motif of larger Mo,S, clusters 
or the bulk MoS; phase. In our calculations, the ground state of Mo2S¢ (C2, 'A1; Fig. 3a) possesses four Mo=S 
double bonds and four Mo-S single bonds, in which each Mo atom has reached its highest oxidation state of 
+6. 

For the neutrals, with the increasing sulfur content, the formal oxidation state of molybdenum increases 
until reaching its highest oxidation state of +6. Meanwhile, the sulfur atoms occupy the terminal sites in Mo2S, 
clusters (n = 4~6) successively. After both Mo atoms get the highest oxidation state, the disulfur units (i.e., S2) 
begin to emerge. Namely, one of the bridging S atoms in Mo2S¢ is replaced by a bridging S2 unit for the Mo2S7 
cluster. Then for the Mo2Ss cluster, two of the terminal S atoms in Mo2S¢ are replaced by the same number of 
terminal S2 ligands. 

Similar evolutionary regularities are also found for the anionic MoS, (n = 4~6) clusters. They preserve 
the structural skeleton of their neutral counterpart. In Mo2S7 cluster, one of the terminal S atoms in Mo2S¢ is 
replaced by a terminal S2 ligand. Subsequently in the Mo2Sg" cluster, two terminal S atoms in Mo2S¢~ are 
replaced by the same number of terminal S2 ligands. Furthermore, based on the relative stability of Mo2S4° 
and Mo2Ss” (Figs. 1 and 5), it seems that the syn configuration is more stable than the anti one. 

4.3 Trend of VDE; as the function of S content in MoS, (n = 4~8) 

Fig. 13 depicts the trend of calculated VDEis as a function of S content in the MoS,” clusters. The trend 

of VDEis: can be qualitatively understood from the frontier MO analysis (Figs. 6~10). The VDEis increases 


nearly linearly as a function of sulfur content (n = 4~6), clearly showing a behavior of sequential sulfidation 
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of Moz dimer. The valence electrons of molybdenum (4d°5s') are sequentially transferred to the added S atoms, 
along with the number of Mo 4d-based orbitals decrease (Figs. 6~8). Interestingly, the VDEis suddenly 
begins to reduce after that. The reason why the VDEj« starts to decrease can also be comprehended from the 
frontier MO analysis as shown in Figs. 9 and 10. When n reaches 7 (i.e., Mo2S7 ), the detachment is origining 
from the fully occupied z* orbitals (71a) of a S2 unit, leading to the decrease of VDEist (from 4.51 to 4.26 eV). 
When n is equal to eight (i.e., Mo2Ss ), we found the VDEsıst corresponds to the detachment from the singly 
occupied z* orbitals (18b1) of two S2 moieties. The further decreasement of VDEis: from 4.26 to 3.58 eV may 
be a consequence of distribution of an unpaired electron over the z* orbitals of two S2 groups. 
4.4 Reduction reaction of Hz on the Mo2S,” (n = 4~8) clusters 

Afanasiev et al. 681 showed that the S2% species located at the edges of fresh MoS; catalysts play a key 
part for the catalytic activity. Recently, Karunadasa and co-workers'**! have reported the synthesis of a side-on 
bound Mo!Y-disulfide complex, which could mimic the MoS; edge sites for the catalytic hydrogen generation. 
It is interesting to note that diverse S2 units are found in the sulfur-rich clusters Mo2S,° (n = 7, 8) in our work. 

The MoS) catalysts are usually used in the gaseous environment of H2/H2S at raised temperature. The 
experiment pointed out that the kinetics of S2 group interaction with hydrogen could be promoted by the 
reactant H? and hindered by the product H2S'©”. Thus, we proposed the reaction (eq. 1) removing a sulfur atom 
from MoS,” (n = 4~8) clusters. The driving force can be estimated by the negative values of Gibbs free 


energy differences (AG) of this reaction, which are obtained by the following equation (eq. 2) and summarized 


in Table 3. 
MoS,” + H2 => Mo2Sy-1” + H2S (1) 
AG = G(Mo2S,-1®) + G(H2S) — G(Mo2S,-””) — G(H2) (2) 


In light of the values of AG, the driving forces of the proposed reactions increase as a function of n. For 
the proposed reaction (Eq. 1; n = 5, 6), the sulfur is removed from the terminal S of MooSs“° and Mo2S6”. 
The corresponding values of AG (Table 3) are positive (Reactant-favored). After both of the Mo atoms reach 
their maximum oxidation state of +6, the S2 units begin to emerge. The energy costs (AG) of sulfur atom losing 
from various S2 units start to be negative (Product-favored). So, for the sulfur-rich clusters Mo2S,,” (n = 7, 8), 
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regardless of the kinetic factor, the proposed reactions (Eq. 2) may be energetically (thermodynamically) 
favored. In the neutral Mo2S7, the bridging disulfide (S2) ligand appears. The Gibbs free energy difference 
(AG) is predicted to be —0.94 kcal-mol!, which should correspond to the removal of S atom from the bridging 
disulfide S2% moiety in Mo2S7. As regards its anion and Mo2Ss¢_", the S atoms from terminal S2% and terminal 
Sx units are removed in the reaction (Eq. 1; n = 7, 8), respectively. The corresponding AG is estimated to be 
-12.51 kcal-mol! (S2* for Mo2S7-) and —14.74 kcal-mol! (S2 for Mo2S¢°), respectively. Generally, the 
driving forces (-AG) that remove a sulfur atom from various S ligands in MoS," (n = 4~8) clusters can be 
sorted in the order: t-S?> < b-S2% < t-S2% < t-S2, which in turn stand for the terminal S atom, bridging 
(2-9!:7') S27, terminal (77) S2% and terminal (4°) S2 unit, respectively. It seems that the S> units are more 
reactive with H2 which would take a S atom away in the form of releasing H2S molecule. Coordination 
unsaturated sites (CUS) can be obtained by removing sulfur atoms from the edges of MoS> catalyst under the 
Hz atmosphere'®’!. The order may provide insight into the the pretreatment of fresh MoS; catalysts under 


hydrogen conditions. 


5 CONCLUSION 


We report a systematical theoretical study on a range of dinuclear metal sulfide clusters: Mo2S, and 
Mo2S, (n = 4~~8). DFT and CCSD(T) calculations were carried out to elucidate the chemical bonding and 


~ clusters. The calculations showed that the sulfur atoms tended 


geometric and electronic properties of Mo2S, 
to occupy the terminal sites of the clusters continuously in the process of sequential sulfidation. After the 
oxidation state of Mo atoms gets the maximum of +6, diverse disulfur ligands emerged in the sulfur-rich 
MoS,” (n = 7, 8) clusters. Additionally, by means of calculating the free energy differences (AG) of the 
reaction (eq. 1), we found that the values of AG were positive for the sulfur-deficient and stoichiometric 
clusters (eq. 2; n = 5, 6), but negative for the sulfur-rich species (Eq. 2; n = 7, 8). The driving forces (AG) for 
the reactions eliminating sulfur from diverse S ligands in MoS, (n = 4~8) clusters followed the general 


order t-S% < b-S2% < t-S2% < t-S2°. This order may provide insight into the pretreatment of fresh MoS 


catalysts under H2 atmosphere. 
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Table 1. Relative Energies of the Low-lying Isomers of the Mo2S,~” (n = 4~8) Clusters 
at the B3LYP Level (AE < 0.40 eV), and Comparisons with Those from 


the CCSD(T) Single-point Calculations at the B3LYP Geometries 


Electronic state B3LYP# > CCSD(T)* © 
Mo2S4 (C3A") 0.00 0.07 
(Cx3A2) 0.04 0.00 
(CB) 0.15 0.38 
(CB) 0.21 0.44 
(C3A") 0.30 0.46 
(Coy!Ay) 0.39 0.22 
Mo2S4- (CyB) 0.00 0.00 
(Con *Be) 0.18 0.17 
(C»2B)) 0.34 0.09 
Mo2Ss; (C/A 0.00 0.00 
(C3A') 0.17 0.40 
MoS5~ (C2A’) 0.00 0.00 
(Cy,4Ap) 0.22 0.41 
Mo2S6 (Cy'A)) 0.00 = 
Mo:5S67 (D2j7Ag) 0.00 Ss 
Mo2S7 (Cy!A) 0.00 0.00 
(C3A") 0.26 0.40 
(C'A) 0.39 0.45 
Mo2S7~ (CPA) 0.00 0.00 
(C2A) 0.39 0.37 
MoS (Cx3By) 0.00 0.00 
(C23 Bu) 0.20 0.23 
Mo2Ss~ (CyB) 0.00 0.00 
(Can Bg) 0.06 0.06 
(C2A') 0.11 0.34 


4 All energies are in eV. °B3LYP/H-BS. ° CCSD(T) single-point calculations at the B3LYP/H-BS geometries 
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Table 2. Vertical Detachment Energies (VDEs) of the Lowest-energy Mo2S,~” (n = 4~~8) 


Clusters and Selected Low-lying Isomer (AE < 0.4 eV) at the B3LYP Level 


MO? VDE» ? Moe VDE»? 
Mo2S4~ 19a;(a) 3.39(Tri) 

(Ca 4B: 0.00 eV) 13b2(a) 3.46(Tri) 
6a2(@) 3.63(Tri) 
18a)(a) 4.51 (Tri) 18a\(f) 3.59(Qui) 
10b,(a) 4.60(Tri) 10b\(B) 4.66(Qui) 
17a\(a) 5.36(Tri) 17a\() 5.18(Qui) 
5ar(a) 5.39(Tri) 12b2(B) 5.23(Qui) 
12b2(a) 5.57(Tri) 5an(B) 5.32(Qui) 
16a\(a) 5.76(Tri) 16a\() 5.65(Qui) 

MooS4~ 9a,(a) 3.31 (Tri) 

(Coy, B, 0.18 eV) 17a,(a) 3.32(Tri) 
15b,(a) 3.81 (Tri) 
16a,(a) 4.34(Tri) 16a,(f) 3.30(Qui) 
Tha) 5.05(Tri) 7b,(B) 4.94(Qui) 
8a,(a) 5.27(Tri) 8a,(B) 4.94(Qui) 
14b,(a) 5.45(Tri) 15a,(B) 5.22(Qui) 
15a,(a) 5.48(Tri) 14b,(B) 5.23(Qui) 
13b,(a) 5.68(Tri) 13b,(B) 5.65(Qui) 

MooS47 19a;(a) 3.32(Sin) 

(Ca 2B, 0.34 eV) 11b,(a) 3.72(Sin) 19a\() 3.25(Tri) 
10b,(a) 4.56(Sin) 10b,() 4.52(Tri) 
12b,(a) 4.90(Sin) 12b2(B) 4.90(Tri) 
5ar(a) 5.07(Sin) 5ar(f) 5.06(Tri) 
18a;(a) 5.23(Sin) 18a\(f) 5.12(Tri) 
17a\(a) 5.43(Sin) 17a\(f) 5.36(Tri) 

Mo2Ss~ 38a'(a) 3.76(Sin) 

(C, ?A' 0.00 eV) 37a'(a) 4.50(Sin) 37a'(p) 4.03(Tri) 
36a'(a) 4.65(Sin) 36a'(B) 4.50(Tri) 
17a"(a) 4.93(Sin) 17a"(B) 4.93(Tri) 
16a"(a) 5.43(Sin) 16a"(B) 5.43(Tri) 
35a'(a) 5.77(Sin) 35a'(B) 5.65(Tri) 
34a'(a) 5.86(Sin) 34a'(p) 5.84(Tri) 
15a"(a) 5.91 (Sin) 15a"(B) 5.88(Tri) 

MooSs~ 5ar(a) 3.81 (Tri) 

(Cx, 4A2 0.22 eV) 25a)(a) 4.11 (Tri) 
13b,(a) 4.62(Tri) 
24a)(a) 4.96(Tri) 13b\() 4.61(Qui) 
13b2(a) 5.22(Tri) 4an(B) 5.21(Qui) 
12b,(a) 5.74(Tri) 13b,(B) 5.21(Qui) 
12b2(a) 5.78(Tri) 12b2(B) 5.73(Qui) 
4az(a) 5.93(Tri) 12b,() 5.77(Qui) 

Mo2S6 15a,(a) 4.51(Sin) 

(Dan 7Ag 0.00 eV) 10b,(a) 4.63(Sin) 10b2,(8) 4.64(Tri) 
Ab3.(a) 4.80(Sin) 4b3,(B) 4.84(Tri) 
8bi g(a) 4.92(Sin) 8bi ef) 4.94(Tri) 
8b1,(a) 5.44(Sin) 8hiAB) 5.43(Tri) 
5bre(@) 5.98(Sin) 5b2,(B) 5.98(Tri) 

Mo2S7- Tla(a) 4.32(Sin) 

(Cı 7A 0.00 eV) 70a(a) 4.75(Sin) 70a(B) 4.26(Tri) 
69a(a) 4.91(Sin) 69a(f) 4.81 (Tri) 
68a(a) 5.00(Sin) 68a(f) 4.96(Tri) 
67a(a) 5.55(Sin) 67a(f) 5.53(Tri) 
66a(a) 5.75(Sin) 66a(f) 5.76(Tri) 

Mo2S7~ 37a(a) 3.53(Sin) 
(Cy 7A 0.39 eV) 34b(a) 4.72(Sin) 34b(B) 4.68(Tri) 
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Mo2Ss~ 
(Coy 2B; 0.00 eV) 


Mo2Ss” 
(Con "By 0.06 eV) 


MoS 
(C,7A'0.11 eV) 


è All energies are in eV. 


36a(a) 
33b(a) 
35a(a) 
34a(a) 
32b(a) 
18b\(@) 
13a2(a) 
27a\(a) 
17b\(@) 
26a;(a@) 
16b\(@) 
21b2(a@) 
15b,(a) 
25a,(a) 
16a,(a@) 
14b,(a) 
23b,(a) 
15a,(a) 
24a,(a) 
52a'(a) 
Sla'(a) 
27a"(a) 
50a'(a) 
26a"(a) 
25a"(a) 
49a'(a) 
48a'(a) 


5.04(Sin) 
5.06(Sin) 
5.45(Sin) 
5.72(Sin) 
5.91(Sin) 
4.42(Sin) 
4.49(Sin) 
4.53(Sin) 
4.80(Sin) 
5.38(Sin) 
5.39(Sin) 
5.76(Sin) 
4,36(Sin) 
4.42(Sin) 
4.61(Sin) 
4.73(Sin) 
5.24(Sin) 
5.38(Sin) 
5.87(Sin) 
4.32(Sin) 
4.63(Sin) 
4.69(Sin) 
4.74(Sin) 
5.37(Sin) 
5.75(Sin) 
5.83(Sin) 
5.98(Sin) 


33b(B) 
36a(f) 
35a(B) 
34a(B) 
32b(f) 


13a2(B) 
27a\(f) 
17bi() 
26a,(f) 
16bi() 
21bs(B) 


16a,(f) 
25a,(f) 
14b,(f) 
23b,(B) 
15a,(B) 
24a,(B) 


51a(B) 
50a'(B) 
27a"(B) 
26a"(B) 
25a"(B) 
49a'(B) 
48a'(B) 


ChinaX ivé (ERAT 


4.97(Tri) 
4.99(Tri) 
5.12(Tri) 
5.68(Tri) 
5.89(Tri) 


3.58(Tri) 
4.52(Tri) 
4.61 (Tri) 
5.23(Tri) 
5.36(Tri) 
5.70(Tri) 


3.73(Tri) 
4.43(Tri) 
4.55(Tri) 
5.07(Tri) 
5.35(Tri) 
5.82(Tri) 


4.24(Tri) 
4.53(Tri) 
4.65(Tri) 
5.35(Tri) 
5.38(Tri) 
5.69(Tri) 
5.90(Tri) 


>The labels “a” and “£” denote the majority and minority spins, whereas Sin, Tri and Qui denote the singlet, triplet and quintet MoS,” (n = 4~8) 


final states upon photodetachment 
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Table 3. Calculated Free Energy Differences (AG) for the 


Proposed Reaction (Eq. 1) at the B3LYP/H-BS Level of Theory 


Losing 


i Losing 
n>n-1 AGne ®? j AGan ©? ; 
sulfur € sulfur € 
534 30.95 t-S?- 39.01 t-S?- 
635 15.45 t-S?- 32.29 t-S?- 
176 0.94 b-S2?- -12.51 t-S?- 
897 -7.04 — 14.74 t-S2- 


a All energies are in kcal-mol!. 
> AGan and AGne represent the AG of anionic and neutral clusters, respectively. 


© The labels “t—S?~’, “b-S?””, “t-S.~” and “t-S»?~” denote the sulfur atom losing from the terminal S, bridging (p2-ņ!:n!) S27, terminal (n?) S77 and 
terminal (n?) S2% upon the reduction reaction (Eq. 1). 
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(g) Mo,S4- Cy, (“By : 0.00 eV) (h) Mo:S7 Cor “By : 0.18 eV) (i) Mo,Sy C», ÊB} : 0.34 eV) 


Fig. 1. Optimized structures for Mo2S4 and Mo2S.4.. The bond lengths are in angstroms (A) 


(a) MoS; C, 1A’: 0.00 eV) (b) Mo,Ss C, ĈA’ : 0.17 eV) 


(c) Mo,S,~ C, ZA’: 0.00 eV) (d) Mo,Sz~ Cy, (A3 : 0.22 eV) 


Fig. 2. Optimized structures for Mo2Ss and Mo2Ss . The bond lengths are in angstroms (A) 
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(a) MoyS¢ Cp, (1A, : 0.00 eV) (b) Mo,S¢- Dan PAg: 0.00 eV) 


Fig. 3. Optimized structures for Mo2S6 and Mo2Sc. The bond lengths are in angstroms (A) 


(e) M02877 Cy ÊA : 0.39 eV) 


Fig. 4. Optimized structures for Mo2S7 and Mo2S7. The bond lengths are in angstroms (A) 
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(e) Mo:S57 G, 7A’: 0.11 eV) 


Fig. 5. Optimized structures for Mo2Ss and Mo2Ss-. The bond lengths are in angstroms (A) 
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(13b,)° Fpsge (19a,)! 
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(19a)! & me Gen 
uv 
e » 
(6a,)! e (6a,)! 
v 
(10b,)}? > (18a,)? 


we 


(18a)? (10b,) 


(a) Mo,S4 Cz, ĈA; 0.04 eV) (b) Mo,S4- Ca, (B; ) 


Fig. 6. (a) Frontier molecular orbitals for the neutral Mo2S4 ground state (Fig. 1b). 
(b) Frontier molecular orbitals for the anionic Mo2S4 ground state (Fig. 1g) 


s " ga , “G66. (38a')! 


of “e = #8 ii 


A 
(36a')? e (36a')? 
© o e 


(a) Mo:S; G, ('A' ) (b) Mo Ss" C, ĈA' ) 


Fig.7. (a) Frontier molecular orbitals for the neutral Mo2Ss ground state (Fig. 2a). 


(b) Frontier molecular orbitals for the anionic Mo2Ss" ground state (Fig. 2c) 


8 8 (25a,)° (15a,)! 
$9 °s 


Se~ Sg- 


(12b)? (dbz)? 


g ē 
(a) Mo:S6 C2, ("Ay ) (b) Mo2S6 Dan Pg ) 
Fig. 8. (a) Frontier molecular orbitals for the neutral Mo2S6 ground state (Fig. 3a). 


(b) Frontier molecular orbitals for the anionic Mo2S6 ground state (Fig. 3b) 
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e 
e 
€ b 


(68a)! 


° (65a)? 


(a) MoS; C) ('A) (b) Mo,S7 Ci ĈA ) 


Fig. 9. (a) Frontier molecular orbitals for the neutral Mo2S7 ground state (Fig. 4a). 
(b) Selected frontier molecular orbitals for the anionic Mo2S7 ground state (Fig. 4d) 


ba © 
* © (27a)? (18b,)! 
= (a d b 
@ ® 
o2. (18b,)? Pa (13a2)? 
òk 
Š % 
E-A m (27a)? 
A T (17b,)! $ (17b,}? 
6° ëg 


(16b;)” } (16b;)” 


(a) Mo Sg Cz, PB ) (b) Mo2Sg~ Cy, ÊB} ) 
Fig. 10. (a) Selected frontier molecular orbitals for the neutral Mo2Ss ground state (Fig. 5a). 


(b) Selected frontier molecular orbitals for the anionic Mo2Ss" ground state (Fig. 5c) 
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AANS 
s S X 
(0.5) (0.5) 
(0.5) 8, » &%.; 
ó » 


(a) MopS C2, CB) 


A iN 

| Js 
4/ X 

(03) (0.3) 

o» 3 8... 
Pd af 

(b) MoS" Cz, 2B,) 


Fig. 11. Valence bond descriptions and numerical electron spin density (in |e|) for the ground state of 


Mo2Ss” clusters. The numerical spin density was shown in parentheses 
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(a) Mo,S,-(C,, *B,) 


(b) Mo,S.-(C, 2A’) 


(c) Mo,S,-(Dy), ?A,) 


(d) Mo,S-(C, 2A) 


(e) MoS- (Cz, *B,) 


04123 4 8 6 
Binding Energy (eV) 


Fig. 12. Simulated photoelectron spectra from the ground states for Mo2S," (n = 4~~8) clusters at the 
B3LYP/H-BS level. The simulations are done by fitting the distribution of calculated VDEs with 


unit-area Gaussian functions of 0.1 eV width 


Fig. 13. Calculated first vertical detachment energies (VDEist) 


of Mo2S,- (n = 4~~8) as a function of S content (n) 
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Probing Diverse Disulfur Ligands in the Mo2S,” (n = 4 ~ 8) Clusters: 


Structural Evolution and Chemical Bonding 


ZHANG Xiao-Fei(#k 5) LIU Xiu-Juan(X7$48) XU Ruo-Nan(4e7% 8) 


WU NI( Ri) HUANG Xin(# WT) WANG Bin(E #6) 


C 


Mox$;- C3, CB) MoS= C, ÊA’) MoS" Dar Ag) MoyS7 CA) Mo$s" C3 ĈB1) 


Density functional theory (DFT) and coupled cluster theory (CCSD(T)) calculations were 
employed to investigate the geometric and electronic structures of a range of dinuclear 
molybdenum sulfide clusters, Mo2S, and MooS, (n = 4-8). After the oxidation state of Mo atoms 
reached the maximum of +6, diverse disulfur ligands emerged in the sulfur-rich MooS,,” (n = 7, 8) 
clusters. The driving forces of removing a sulfur atom from different S ligands in Mo2S,,”° (n = 
4~8) clusters, especially from those disulfur units, were evaluated. Vertical detachment energies 
(VDEs) were predicted, and then the photoelectron spectra (PES) were simulated. Molecular 


orbital and spin density values were analyzed to elucidate the chemical bonding and the 


evolutionary behavior in the dinuclear molybdenum sulfide clusters. 
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